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COMPRESSIBLE  TURBULENT  BOUNDARY  LAYER 


Bao  Hanling 


Chinese  Aerodynamic  Research  and  Development  Center 


^  This  paper  presents  an  improved  en¬ 
trainment  method  for  the  two*dimensional 
compressible  turbulent  boundary  layer. 

This  method  is  simple  with  short  computer 
time  and  can  be  used  for  the  development  of 
calculations  for  the  turbulent  boundary 
layer  behind  a  shock  wave.  Calculated 
results  from  six  examples  indicate  that  this 
method  has  excellent  numerical  accuracy. 

When  compared  with  P.  D.  Smith's  method,  the  ___ 
results  show  a  marked  improvement.  1 C..-  -  . 
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1.  Preface 


The  calculations  of  turbulent  boundary  layers  at  aircraft  wings 
for  the  modification  of  flow  field  of  transonic  nonviscous  flow  is 
indispensible .  Although  many  differential  methods  have  been  devel¬ 
oped  abroad,  certain  practical  transonic  computer  programs  (e.g.,  the 
three  NASA  full  potential  programs:  FL-023,  FL-027  and  FL-030)  still 
utilize  the  integral  method  as  the  counterpart.  This  is  because  the 
integral  method  is  cost  effective,  convenient  to  use,  and,  at  least 
at  present,  not  inferior  to  the  differential  method  in  accuracy. 
Currently,  the  entrainment  method  appears  to  be  the  only  integral 
method  that  can  be  extended  into  the  application  of  three-dimensional 
aircraft  wings.  This  method  was  developed  from  the  basis  of  M.  R. 
Head's  mechanical  model  by  P.  D.  Smith^^.  Based  on  our  calculating 
experience,  the  accuracy  for  certain  examples  using  the  mechanical 
model  of  P.  D.  Smith  is  not,  even  for  two-dimensional  condition, 
satisfactory.  This  paper  attempts  to  improve  this  method.  Since 
examples  used  for  evaluating  accuracy  of  the  three-dimensional  turbu¬ 
lent  boundary  layer  are  extremely  limited,  this  paper  only  validates 
the  effects  of  the  improved  method  for  two-dimensional  conditions. 

r  21 

M.  R.  Head  and  V.  C.  Patel  were  the  first  to  present1  J  the 

method  of  taking  the  effects  of  the  non-equilibrium  boundary  layer 

into  consideration  to  improve  the  mechanical  model  of  the  entrainment 

method.  Their  method  only  applied  to  a  two-dimensional  incompressible 

flow.  In  1973,  J.  E.  Green,  et  al  developed  the  lay  entrainment 
r  3 1 

method1  i.e.,  to  use  one  differential  equation  to  determine  the 
entrainment  coefficient  Cg,  and  thereby  taking  the  upstream  effects 
during  entrainment  into  account.  This  method  has  been  extensively 
applied  abroad,  yet  so  far  no  reports  have  been  made  public  about 
extending  it  into  the  calculations  for  three-dimensional  boundary 
layers.  This  paper  adopts  another  approach  and  attempts  to  extend 
the  improved  method  of  Head  and  Patel  to  compressible  boundary  layer 
conditions.  A  method  which  takes  the  lag  effects  of  the  pressure 


gradient  into  consideration  to  improve  the  calculations  for  the 
entrainment  coefficient  under  a  strong  adverse  pressure  gradient  is 
also  presented  in  order  to  calculate  the  development  of  the  turbulent 
boundary  layer  behind  a  shock  wave.  The  results  of  six  examples 
indicate  that  the  improvement  effects  of  this  method  are  not  inferior 
to  Green's  method  and  that  the  calculation  accuracy  is  also  not 
inferior  to  a  complicated  differential  method. 

2.  General  Description  of  Method 


The  two-dimensional  boundary  layer  entrainment  method  must  solve 
the  momentum  integral  equation  and  entrainment  equation  simultaneously, 
and  they  are 

2  +W“‘W»)"C|/ 2  (i) 

( 2) 

where  0^  is  the  momentum  thickness;  C^,  is  the  entrainment  coefficient; 
Cj  is  the  wall  surface  friction  coefficient;  H  and  are  the  shape 
factor  and  entrainment  shape  factor;  L  is  the  pressure  gradient  param¬ 
eter:  — ^jf**"*  ant*  ^e  are  t*ie  velocity  outside  boundary  layer 

i  d  Mach  Number,  respectively.  Equations  (1)  and  (2)  can  be  closed 
by  substituting  several  emphirical  and  semi-emphirical  formulae  re¬ 
lated  to  the  above  parameters.  The  calculated  results  of  the  boundary 
layer  are  greatly  affected  by  how  C£  is  determined.  In  P.  D.  Smith's 
method,  is  completely  dependent  upon  the  local  boundary  layer's 
total  parameters,  and  thus  the  history  of  upstream  effects  of  tur¬ 
bulent  flow  cannot  be  considered.  Based  on  the  thoughts  of  Head  and 
Patel,  this  paper  expresses  CE  as 

C.-C.^F(f)  (3) 

where  represents  the  value  of  the  boundary  layer  under  local 

equilibrium  conditions;  F(r)  represents  the  non-equilibrium  factor 
which  shows  the  deviation  situation  of  C-,  under  equilibrium  conditions; 

F  r  o  1 

r  is  the  characteristic  parameter.  Head  and  Patel  pointed  out1  J 


where 


r<!  F(r)-(5 -Ar  )/(S-*r) 

'>1  F(  r  )-  1/(2  r  -  l  )  (5) 

Morkovin^^  assumed  that  turbulent  flow  structure  is  not  affected  by 
the  Mach  Number  as  long  as  the  Mach  Number  based  on  typical  pulse 
velocity  is  rather  small.  Therefore,  we  can  regard  F(r)  as  unchanged 
for  compressible  flow  and  only  rewrite  r  as 


1  d  (P«*«8u)  /f  1 

dx  /  Ip««,  dx  J. 


(6) 


(here  Pe  is  the  flow  density  outside  the  boundary  layer), 
equilibrium  formula  can  be  rewritten  as 

<?-C4  82S5(*  +1.0717),/*+l.»«iX  1  +0.04Me,),', 

where 


The 


G-(h~  I  )/  «/ 

•--(tH/C,)  a*n)« 


(7) 

(8) 

(9) 


Here  H  is  the  conversion  shape  factor.  Utilizing  equations  (1),  (2) 
and  (6)'*''  (9)  as  well  as  the  formulae^*]  for  H,  and  R,  it  is  not 
difficult  to  solve  for  r  and  CEe(^,  and  therefore  CE  from  equations 
( 3 )  and  ( 5 ) . 


When  a  shock  wave  is  present  on  the  wing  surface,  the  effects 
of  a  pressure  jump  across  the  shock  wave  will  sustain  far  into  regions 
behind  the  wave.  Therefore,  in  order  to  take  these  effects  into 
account,  it  is  appropriate  to  substitute  L  in  equation  (6)  with  the 
lag-pressure  gradient  parameter  L^.  Inspired  by  Reference  [6]  and 
applying  the  functional  form  of  lag  model's  weighting  function  there¬ 
in,  the  expression  for  we  substituted  was 

L> "4“  *  -jr.mti./fW*, 

where 

0  " 51  e,pC_(  *  -*>y/MlH'-)dx, 


(10) 

(11) 


xn  represents  the  starting  point  of  the  shock  wave;  \  is  a  constant 


associated  with  the  pressure  jump  value  and  based  upon  experience 


it  can  be  expressed  as 


X  -185AC, 


(12) 


3.  Calculated  Results 

Six  examples  of  a  profile  with  a  compressible  turbulent  boundary 
layer  Mach  Number  range  of  subsonic  and  transonic  flows  were  calculated 
using  the  present  method.  The  pressure  distribution  used  in  the  cal¬ 
culations  were  measured  in  a  wind  tunnel. 

Examples  1  and  2  are  the  boundary  layers  for  upper  and  lower 
surfaces  of  RAE  2814  profile.  In  Figs.  1  and  2  the  calculated  results 
are  compared  with  experimental  results  (M.aar0.725,  R*— 1.5X10*).  It  can 
be  observed  from  the  figures  that  the  results  of  P.  D.  Smith's  method 
and  those  of  the  experiments  for  the  upper  surface  of  the  profile 
differ  little.  They,  however,  differ  very  much  for  the  lower  surface 
of  the  profile.  Results  of  this  method  show  good  agreement  with  ex¬ 
perimental  data  for  both  upper  and  lower  surfaces.  For  comparison 
purposes,  calculated  results  of  Green's  method  are  also  given  in 
Figs.  1  and  2.  It  can  be  observed  that,  for  these  two  examples,  the 
improvement  effects  of  Green's  method  is  inferior  to  the  present 
method  and  its  initial  values  for  the  upper  surface  of  the  profile 
show  greater  deviation  from  the  experimental  points.  Moreover,  Green's 
method  is  more  complicated  than  the  present  method  and  consequently, 
computation  time  is  longer.  The  experimetal  results  for  this  example 
are  from  Reference  [3]. 
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Figure  3  shows  the  two  sets  of  calculated  results  for  the  upper 
and  lower  surfaces  of  RAE  2822  profile  using  the  present  method  and 
P.  D.  Smith’s  method:  (a)  Moo=0.676,  Re=5.7xl06,  oL  2.18°;  (b)  Mw  = 
=0.73,  Re=6.5xl0^,  0(=3.19°.  For  condition  (a),  no  shock  wave  is 
present  on  the  wing  surface  and  data  from  the  present  method,  P.  D. 
Smith's  method  and  the  experiments  are  in  good  agreement.  For  con¬ 
dition  (b),  a  shock  wave  is  present  on  the  wing  surface.  The  cal¬ 
culated  results  of  and  H  using  Smith's  method  show  marked  devia¬ 
tion  in  the  region  immediately  following  the  shock  wave,  whereas  those 
using  the  present  method  are  basically  consistant  with  the  experimen¬ 
tal  results,  as  each  H  displays  a  peak  value  following  the  shock  wave. 
The  experimental  results  for  this  example  are  from  Reference  [4]. 


Figure  4  is  the  calculated  results  for  the  low  speed  boundary 
layer  of  the  NACA-65( 216)-222  profile.  The  attack  angle  #=10.1°. 

The  calculated  results  using  the  present  method  and  the  Reynold's 
stress  equation  as  well  as  the  experimental  data  are  compared.  Even 
though  the  Reynold's  stress  equation  is  more  stringent  in  theory  than 
the  present  method,  its  computation  time  is  much  longer  and  the 
calculated  results  are  not  necessarily  better  than  those  from  the 
present  method.  Except  for  points  near  the  separation  point,  the 
results  of  the  present  method  are  very  consistent  with  the  experimental 


data.  Both  theoretical  methods  predicted  the  separation  trend  after 
-——>0.55  . 
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Fig.  4.  Calculated  results  for  upper  surface  of  NACA-65 ( 216 ) -222 
profile  0<=10.1o,  Re  =  2.6xl06 

Key:  (1)  Reynold's  stress  equation  method;  (2)  Experiment; 

(3)  Present  method. 


Figure  5  gives  the  calculated  results  for  the  upper  surface  of 
the  super  critical  profile  DSMA-523.  The  calculation  conditions  are: 
Moo=0.83,  Re=3xl0^,  0^=0. 54.  A  shock  wave  is  also  present  on  the 
wing  surface.  The  calculated  results  of  this  example  are  not  as  good 


as  those  of  the  previous  five  examples,  and  this  is  related  to  the 
characteristics  of  the  super  critical  profile.  The  results  of  T. 
Cebeci's  differential  method  are  also  given  in  the  figure.  Though 
this  method  consumes  for  more  computation  time  than  the  present  method 
its  calculated  results  are  not  necessarily  more  accurate  than  those 
of  the  present  method.  Especially  for  the  shape  factor  H  this  method 
fails  to  show  the  peaks  which  appeared  following  the  shock  wave  and 
the  rising  trend  near  the  tail  end. 


Fig.  5.  Calculated  results  for  the  upper  surface  of  DSMA-523  profile. 
M oo=0 . 83 ,  Re=3xl06,  Cy=0.54 

Key:  (1)  Present  method;  (2)  Cebeci's  method;  (3)  Experiment. 

4.  Conclusions 


This  paper  presents  an  improved  entrainment  method  for  calculating 
a  two-dimensional  compressible  turbulent  boundary  layer.  When  there 
is  a  strong  adverse  pressure  gradient  on  the  wing  surface,  especially 
when  a  shockwave  is  present,  the  lag  effects  of  the  pressure  gradient 
near  the  shock  wave  are  taking  into  account  and  rather  satisfactory 
calculation  results  can  still  be  obtained.  Calculated  results  from 


examples  of  six  profile  boundary  layers  show  that  the  accuracy  of  the 
present  method  compares  favorably  with  those  of  some  of  the  outstand¬ 
ing  calculation  methods  for  a  turbulent  boundary  layer.  Also,  the 
present  method  is  simple  and  fast  and  easier  to  extend  to  three- 
dimensional  cases.  The  mechanical  model  of  the  present  method  has 
currently  been  introduced  into  the  calculation  of  a  three-dimensional 
turbulent  boundary  layer  on  aircraft  wings. 

Special  thanks  to  Professor  Bian  Yinguei  of  the  Institute  of 
Mechanics,  Academia  Sinica,  who  reviewed  and  advised  on  this  paper. 
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(A  new  improved  entrainment  method  for  calculation  of  compressible 
♦urbulqnt  boundary  layer  is  suszcsted,  in  which  the  effect  of  non-equilib- 
r  u  j  turbulent  boundary  layer  on  determining  entrainment  coefficient  is 
taken  into  account,  and  a  kind  of  lag-pressure  gradient  parameter  is  intro¬ 
duced  into  calculation  of  flow  field  with  strong  adverse  pressure  gradient, 
such  as  that  near  shock  waves.  For  .evaluation  ind  validation  of  the  me¬ 
thod  six  numerical  examples  are  presented,  and  the  results  show  good 
agreement  with  experimental  data.  Moreover,  in  comparison  with  ether 
threoretical  methods  the  present  method  is  aa  good  as  J.  E.'  Green’s 
method,  T.  Cebeci's  method  and  Reynold's  stress  equation  method,  but 
more  efficient  in  computation,  aad  suoerior  to  I*.  D.  Smith's  method  in 
approaching  to  the  experimental  data./The  method  is  also  extended  to 
threedimeasioaal  cases. 
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